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ABSTRACT 

The chemical composition of cosmic ray nuclei with 3KZ<28 between a few 
hundred MeV/nuc and a few hundred GeV/nuc are compared with a consistent set 
of propagation calculations. These include the effects of spallation (energy- 
dependent cross sections are used), escape, ionization loss in the inter- 
stellar medium and deceleration in the solar cavity. The amount of matter 
traversed by cosmic rays is found to be about 7 g/cm a , independent of energy 

between 100 MeV/nuc and 2 GeV/nuc. Above 2 GeV/nuc, the escape length varies 
0 4tfc0 • 1 

as E“ * ' ’ . In addition, a procedure has been developed to measure the 
shape of the cosmic ray pathlength distribution. 

Utilizing the ratio of Fe-secondaries to Fo in the cosmic rays, -presently 
available data are found to be consistent with an exponential distribution, 
and to eliminate models in which the path length distribution is severely 
truncated. To tie down the shape of the distribution more precisely, new 
measurements of the cosmic ray composition presently becoming available from 
experiments on the HEAO-3 satellite will have to be coupled with improved 
measurements of the energy dependence of partial and total cross sections. 

Subject Headings ; cosmic rays: abundances — interstellar: matter 
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Introduction 

Hie chemical composition of cosmic rays reflects information not only 
on the mean amount of matter traversed by the particles between their sources 
and the solar system, but also on the shape of the distribution of matter as 
sampled by the particles. As is well known, cosmic ray nuclei from sources 
(primaries) spall on traversing the matter in interstellar space, producing 
secondary nuclides. By studying the relative composition of heavier 
secondary nuclides which are rare or absent at the source, the shape (or 
moments) of the path length distribution can be studied. This shape reflects 
the location (either spatially or temporally) o ‘ sources relative to the 
observer (Owens, 1976), and hence Is important ,n understanding their origin. 

In this paper We present calculations of t le effects of the distribution 
of path lengths through which the cosmic rays p 'opagate on the detailed 
chemical composition and compare these calculat .ons with data from a variety 
of experiments ranging from 100 MeV/nuc to 100 5eV/nuc, 

It has been generally believed that the co position was beat explained 
by a path length distribution (PLD) with an abs mce of short path lengths. 

This was thought necessary to explain the large number of nuclei observed in 
the sub-iron region (Shapiro and Silberberg, 19/0; Shapiro, Silberberg and 
Tsao, 1973; Gatcia-Munos and Simpson, 1970; Fic itel and Reames, 1968). Since 
that time new measurements have been made of so lie of the cross sections 
involved. Hie mean matter traversed has been f mnd to vary with energy, and 
more accurate charge composition measurements hive been made. 

There is a sound theoretical basis for a P JD in which short path lengths 
are absent. Lesniak and Webber (1979) have derived the PLD shapes expected 
in some of these models. In one such model somi of the matter through which 
cosmic rays pass is located in the cosmic ray source. Cowsik and Wilson 
(1973, 1975) first Invoked a ’’double or nested leaky box” model to explain the 


// 
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early observations of energy dependent changes In the abundance ratios of 
primary cosmic rays. Simon (1977) Introduced a convenient parameterization 
of the PLD resulting from these models in terms of the convolution of two 
exponential distributions. We utilize this formalism to Investigate the 
shape of the PLD using the cosmic ray* data. 

In this paper, the available data will be compared with a self consistent 
set of calculations. By omparing the data from different groups with the 
same model, any differences which might exist between calculations are 
eliminated in interpreting the data. 

We will assume that acceleration of cosmic rays during their propagation 
in the interstellar medium is negligible and that the cosmic rays are effec- 
tively Injected into the galaxy with a power law spectrum common to all 

species. This would apply if cosmic rays are accelerated in localized 

accelerators and then injected into the interstellar medium where the probability 
of subsequent acceleration is small. This could be localized shocks (Blandfqrd 
and Ostrlker, 1978; Jokipli and Higdon, 1979), the acceleration of a thermal 
plasma to cosmic ray energies (Eichler, 1980) or acceleration in supernova 
remnants (Scott and Chevalier, 1975). 

We begin with the simple leaky box model (Cowsik et al. 1967), a model 

of cosmic ray propagation in which the PLD is exponential. Using the observed 

relative abundance of the light secondary elements Be and B and the primary 

elements C and 0, we can determine a mean escape length X (E) as a function 

e 

of energy. The general featu -es found by Ormes and Freier (1978) are present: 
a decreasing escape length at high energies (> 2 GeV/nuc), and a constant escape 
length at lower energies. Jones (1979) has shown that the general shape of 

•I ■ 

X p (E) can be explained by a galactic wind model (Joklpil, 1976) in which escape 
of particles by diffusion dom: nates at high energies while escape by convection 
dominates at low energy. 
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The effect of varying the shape of the PLD on the secondary nuclei in 
the range 2KZ<25 just below the abundant primary iron nuclei is studied since 
the abundance of these elements is more sensitive to the shape of the PhD than 
the light secondary elements. Such studies are dependent upon a knowledge of 
the partial cross sections for producing each secondary isotope from every 
possible primary Isotope. Insufficient measurements of the energy dependence 
of these important cross sections have been made as noted by Raisbeck (1979) 
and care must therefore be taken when using these data. We use here the semi- 
empirical formulations of Sllberberg and Tsao (1973, 1973a, 1977, 1977a, 1977b 
and Tsao and Sllberberg, 1979) including all their corrections based upon the 
recent measurements with the Bevalac. The various .versions of these formu- 
lations have been cross checked with Sllberberg and Tsao under a community- 
wide effort initiated by P. S. Freler. We shall discuss the limitations on 
the interpretations due to uncertainties in the cross sections. 


Method of Calculation 

In the equilibrium steady-state homogeneous model, the intensity 

J^(E) of particles of type 1 at energy E is given by 

J,(E) J . (E) J.(E) J.(E) 

q ( E) - -i 1 1 + E { — U 

V E > X dec< E > X lnt (E) “i 1 X spall< E) 


J k (E) 


> - ^ [w i (E)J i (E)] - 0 


( 1 ) 


1 i 1 

where is the source term; X^. ^dec ant * ^int are t * le mean P at ^ lengths 

• ki ki 

for escape, decay and interaction of nuclei of species i; and \j ec 

are the mean path length for spallation and decay of nuclei of species k 


into species 1; -w^(E) is the mean energy loss per unit path length. Note 


that X 


<m> 


int a. 

(sp* 11 ) {spaU) 


and X dec ■ yT dec $cn<m>, where n and <m> are the mean 


number density and mass of interstellar nuclei encountered by the cosmic 

-3 

rays, taken to be 0.3 cm and 1.3 amu, respectively. 
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In actual calculations, solutions J^ B (E,x) to the slab model equations, 
^ J s (E,x) J, s (E,x) J S (E,x) J k S (E,x) J k S (E,x) 


+ ' 5 e [ w i< E > J i 3(E ' x):! * 0 • 


with initial conditions, J^ s (E,0) - (^(E) , may be obtained numerically; 

j s (E,x) - j ^757“ r>1 r f- J 

1 A 1 i (E> hit™ hL™ 


J. v W k (®«) s 
+ i v. J, S (E, ,x-Ax) 

k>i w, (E) k v k’ ' 


Call® X d£® 


Here E is the energy of a nucleus of species i at pathlength (x-Ax) such 
that its energy at pathlength x is E. E^ is obtained from range-energy 
tables or, for small relative losses in energy, E^ « E-w^CEJAx* 

The solutions J^E) to the full equation may then be obtained from 
J S (E, X ) by weighting with the pathlength distributions P^(x»E) (e.g* 

Leznlak, 1979) ; 

00 

J^(E) * / J i a (E,x)P jL <X,E)dx (A) 

For an exponential PLi), P. (x,E) 53 expC-x/\^(E)3 . For cross checking purposes, 
when using an exponential pathlength distribution, we can Include an additional 
term -Ax/\ e (E) within the first curly bracket of equation (3) and use P^(x,E) ■« 1 
in equation (4) * 

For the spallation cross sections, we have used the Sllberberg 

and Tsao semi-empirical formulae and taken an interstellar medium 

consisting ot 90 percent hydrogen and 10 percent helium by number. Some of 

the Important cross sections at just a few of the energies used in our calcula- 

12 

tlon for the spallation of C and Fe are given in Tables 1 and 2. These 
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cross sections include thr> effects of the decay of isotopes with half lives 
less than 10 s and should thus be considered the "cosmic ray cross sections « M 
In each case the cross section given la in mb pur average interstellar nucleus 
and for Table 2 the cross sections have also been weighted with the expected 
abundances of Fe Isotopes in the cosmic radiation* These tables Illustrate 
the importance of using the observed energy dependence of the crosB sections; 
e*g*, the cross section for production of Gr by spallation of Fe varies by 
a factor of 4 between 100 MeV/nuc and high energies. Any error in the 
spallation cross sections Will be directly reflected as an error in that 

predicted secondary- to-primary ratio. 

* 

Ionization losses are also very important in determing cosmic ray spectra 
at low energies and cannot be Ignored. The effects of neglecting the energy 
dependence of the spallation cross sections and of neglecting Ionization losses 
on the ratios of Fe-secondaries to Fe and (Be+B) to C are illustrated in Figure 1 
(for the model we adopt as our "standard" to be described later). 

The total cross sections adopted also affect the predicted composition, 
particularly for heavy nudlei (e.g„, sub-Fe/Fe ratio) where particle losses 
are dominated by nuclear interaction rather than escape. We use the energy 
independent fits to the da\!a obtained by Hagen (1976) : 

a int " 3.0 Tr <l‘29) i CA T 1/3 +A B 1/3 - 1.189 exp{-0. 05446 min(A T , A fi )}] 2 mb (5) 

where A^, and are the mass numbers of the target nuclei (interstellar matter) 
and the beam (cosmic tay nuclei). An energy dependent total cross section scaled 
from the (p, p) cross sectior (Karol, 1975) has been used by Garcia-Munoz et al. 

(1979) in their analysis* While perhaps theoretically justified, the present 

* 

data are unable to determine the energy dependence with any certainty and new 
measurements are urgently needed* The values used here (e.g., 716 mb for the 

, total inelastic cross section) may be somewhat low. Values in the literature 
range from 721 A 7 mb (Roberts et al. 1979) to 779 * 27 mb (MacFall et al., 

1979) for measurements at energies above 100 GeV, but the inelastic cross 
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section rises with energy, The former cross section Is an n-Fe measurement and the 
latter a p-Fe measurement. What Is needed In our calculations Is the mass changing 
cross section (taking account of Fe •* Fe + n) and the only measurement is 750 & 50 mb 
(Westfall et al,, 1979) * Our adoption of the Karol (1975) energy dependent form 
would reduce our predicted Sub-Fe to Fe ratio by about 10 percent above 400 HoV for 
the same grammage. 

In our calculations we use the elemental source composition derived by Leznlak 

and Webber (1978) In conjunction with the Isotopic abundances given by Cameron (1980). 

These are given in Table 3. Elemental abundances for Cl-K and Sc-Mn are alBo from 

Cameron (.1980) but normalized to the upper limits given by Luznink and Webber (1978) 

«2 3 

for these groups. All injection spectra are taken to be proportional to W * where 
W la total energy per nucleon (Ormes and Freler, 1978). Hie abundance ratios 
calculated depend on the Injection spectra and primary composition adopted. 

We have Included the effects of solar modulation on the observed composition 
by using the "Force Field" solution (Gleeson and Axford, 1968) , The modulated 
spectra J^ m (E) are obtained from the Interstellar spectra J^(E) at a given 
kinetic energy per nucleon E by 



E 2 +2E 

(EH) 2 +2(E+*) 


J x (E+*) 


( 6 ) 


Ze 

where $ *• is the mean energy loss per nucleon in the heliosphere, We have 
usually used a deceleration parameter <{> " 200 MV as being typical of values 
derived from the data, although the effect of varying t|> is discussed later. 

We also consider the effect on composition of propagation through the 
atmosphere. For this we solve the slab model equations (equation 3) with initial 
conditions J i s (E,0) « J^CE). For the case of an experiment at a mean slant depth 
x^ in the atmosphere we take a pathlength distribution P(x) * 6 (x-x^) . In 
reality, this will be broadened by the finite acceptance solid angles of detectors; 
but, since this varies from experiment to experiment, we do not consider this 
effect here. The amount of atmospheric correction necessary can be seen in Figure 2. 
Here we present results of our calculations of the ratios (Be+B)/C and (2KZ<25)/Fe 
in interstellar space, abovn the atmosphere (effects of solar modulation with 


8 


<(> » 20G MV) and ut various depths In the atmosphere appropriate to balloon-borne 
experiments, these calculations are for our standard model (to be described in 
the next section). As can be seen, features in the ratios below about 1 GeV/nuc 
will be washed out in a few g/cra of atmosphere because of ionization losses. 
Interpretation of Be, B, C Ccmpositlon 

i 

The role of the ratio cf (Be+B)/C is to determine the mean escape 

length as outlined in the introduction. Here we present the results of our 

calculations of the ratio (Be+B)/C ub a function of energy for an exponential 

pathlength distribution,. The results are shown iti Figure 3 for various mean 

escape lengths,. Here we have used a solar modulation deceleration parameter of 

<j> * 200 MV. We also show a survey of the data on this ratio corrected to 

the top of the atmosphere (in most cases by the original authors). These 

predictions have been used tc derive a mean escape length for each data 

point and Figure A shews the variation of mean escape length with energy 

obtained in this way. Above 2 GeV/nuc the data are generally interpreted 

«s suggesting that the mean escape length is proportional to E 

Below 2 GeV/nuc we find that the mean pathlength is essentially Independent 

2 

of energy at g/cm . This conclusion is based on the data at the lowest 
energies considered (M.Q0 Me v /nuc) and is thus sensitive to the amount of 
solar modulation adopted (<J)-200 MV) and to the shape of the injection spectrum. 

» *s 

The data at ~1 GeV/nuc, where solar modulation is not **qo Important, fix the 

2 

escape length at about 7 g/cm ; but at energies around 100- MeV/nuc an escape 

2 2 
, length as high as 9 g/cm (no solar modulation) or as low as 5 g/cm 

k 

(<j> « 400 MV) is also possible 1 . For our "stande rd model” we therefore 
adopt 


V E > 


7 g/cm 2 

nt E n - 0.4 ,2 

2GeV/nuc * 


E<2 GeV/nuc 
E>2 GeV/nuc 


(7) 


for all species with Z>2. Predictions of the ,*atio (Be+B)/C based on this 
standard model are shown by the dashed curve o E figure 3. 
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The energy dependence of the mean escape length, above Is not expected 


to be valid at the highest energies. For example, the mechanism responsible 
for this energy dependence woild almost certainly break down at energies for 
which equation (7) would give a mean esffjftpe length less than the thickness of 

i 

the galactic disk. Some information about the high energy behaviour may be 

p 

» 

obtained from measurements of the anisotropy of cosmic rays. Provided the 


* streaming of cosmic rays Is n>t much In excess of that necessary to transport 
the cosmic ray nuclei from the sources to the escape boundary, then the magnitude 
of the expected anisotropy Is given by (Klraly et al. , 1979): 


t con 


where t con Is the confinement time of cosmic rays In the galaxy and l* DE is the 

"direct exit" time, i.e., the escape time of neutral cosmic rays. Since the 

confinement time is proportional to escape length, the anistropy would be 

inversely proportional to escape length. The treasured anisotropy is almost 

o 

constant at about 0.05 percent over the energy range 5 x 10 GeV/nucleus to 
5 x 10^ GeV/nucleus (see the review by Elliot, 1979), and this has been used 
by Kiraly et al. (1979a) to Infer that the mear escape length is essentially 
constant over this energy range. On the other hand, the absence of any 
flattening in the cosmic ray proton spectrum up to 2 TeV/nuc suggests that 
the energy dependent escape irechanisms continues to at least that energy. 

A flattening in the all particle spectrum between 10 and 100 TeV is seen 
in air shower spectra. Perhaps this flattening is associated with a change 
in the escape mechanism. 


The abundances of the (dements Li through N1 predicted at various 
energies for the standard model are given alongside the source abundances 
in table 3. Using this model we can look at other abundance ratios to see 
whether they are consistent with the experimental results. In figure 5 
we reproduce the survey of data on the ratio (175.2^25) /(26^Z< 28) from Simon 



et al, (1980). Tills Is compare 1 with the ratio predicted for the standard 
model In the figure. The energ; dependence of the predicted ratio Is in good 
agreement with the observations , bearing in mind the statistical uncertainties. 
In particular , we are able to reproduce the observed peak in the abundance 

i * 

ratio at ** 1 GeV/nuc. 

From the variation of the ratio of two "primary" components with energy 
we can look for variations in the slope of injection spectra between the 
different elements. A survey of vie^aureroents of the Fe/0 ratio is given in 
Figure 6, and it is seen that, despite the large amount of scatter between 

* * i 

the data, this ratio increases from ~ 0,8 at 100 MeV/nuc to ^ 1,5 at 100 

GeV/nuc. The predicted ratio also increases with energy due to t\v. energy 

dependence of the relative importance of escape and interaction losses, 

especially for Fe nuclei. At \ery high energies, since losses are dominated 

by escape (X e <<X lnt ) for both species, the Fe/O ratio approaches the source 

ratio. Below about 1 GeV/nuc, the ratio is also reduced by ionization losses. 

These effects ate further enhanced for observations made from balloons in the 

atmosphere. The predicted ratio is shown in Figure 6 for depths of 3, 6 and 
2 

9 g/cm in the atmosphere. The scatter in the data is quite large and could 
arise partly from the use of different atmospheric correction factors by the 
various authors, Within the uncertainties in the data we do not detect any 
significant difference between the injection spectra of Fe and 0; in this 
energy range the observed flattening can be explained solely by propagation 
effects as noted by Cowsik and Wilson (1973). 
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^'Two-zone 11 or "Neatod Leal^ v, Box 1 *' Mbdela 

The composition of heavy nuclei is particularly oensltive to the 

shape of the pathiength distribution because the interaction mean free’ paths 

Pe 2 

of heavy nuclei are much less than the mean escape length <Q.g* ^int 82,6 S /cra >* 

Thus if the pathiength distribution were truncated for small pathlengtha, 

then nuclear interactions would have a greater effect in determining the 1 * * 

observed composition, resulting in a significantly higher ratio of secondary 

nuclei to primary nuclei. A similar effect, although not so pronounced, 

results for the composition of light and medium nuclei, increasing,, e.g. , the 

(Be+B) /C ratio observed, and this must be taken into account when deriving 

the mean escape length from the data for such models, 

Cowslk and Wilson (1973, 1975) have shown that "two-zone” or "nested 

leaky ’box” models result in a PLD which is the convolution of two exponential 

distributions. We use the formalism Introduced by Simon (1977) to give 

physical meaning to a parameter which describes the shape of the PLD. The 

distribution is uniquely defined by two parameters: the- mean pathiength in 

the source region X g , and the mean pathiength in the galaxy X^ (the mean 

escape length is X g » X g +X b , ind the ratio* 5 X R /X e determines the shape of 

the distribution). Pathlengtn distributions for X g /X fi ranging from 0.0 

(pure exponential PLD) to 0.5 (for X /X, values greater than 0.5, the shape 

of the distribution is the sane as for 1.0-X /X,) are shown in Figure 7 for 

s a 

2 

X *** 5 g/cm . We use the ratio X /X tt here as e convenient parameterization 

6 S vJ 

of the shape of the PLD. Once the. shape has been determined experimentally, 
it can be compared with different cosmic ray storage models (Leznlak and Webber, 
1979) or source distributions (Owens, 1976). 

The variation of the \B/C and (Be+B)/C ratios with mean escape length X £ 
is given in Figure 8 for various values of X g /X e . The calculations for this 
figure were done using the asymptotic cross sections and correspond to energies 
above 4 GeV/nuc. The cross sections are however almost Independent of energy 
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above 1 GoV, This Is fortunate since the (Bc+B)/C ratio is best known at around 

1-2 GeV/nuc and we ahall uae this to fix X for different values of X /X. 

e a e 

From figure 8 we see that X ft could vary from 5.2 g/cm 2 <X 8 /X- e - 0.5) to 
2 

7.0 g/cra (X /X “ 0), depending on the shape of the PLD, and give the 
s e 

observed (Be+B)/C ratio between 1 and 2 GeV/nuc. 

We have calculated the ratio Of Iron secondaries (21<Z<25) to Fe and 
(Be+B)/C as a function of energy for two different cases; the standard 
model (X /X - 0; \ - 7 g/cm 2 below 2 GeV/nuc, 7(E/2)”®‘^ g/cm 2 above 

0 U 6 

2 

2 GeV/nuc) and for a eeverly truncated PLD (X /X « 0.33; X * 5.3 g/cm 
below 2 GeV/nUc, 5.3(E/2)~®*^ g/cm 2 above 2 GeV/nuc). The predictions ore 
shown in Figure 9. In the figure, we also show the effect of varying the 
amount of solar modulation and give results for <}>' * 0, 200 MV and 400 MV. 

The predictions for the (21<Z<25)/Fe ratio may be compared with a survey of 
the experimental r« *ults shown in the figure. From this comparison, it 
ap that there is no need to invoke a truncated pathlength distribution 


\ 


& 


to explain the observed ratio. The data are consistent with the predictions fpp the 
exponential PLD; however, because of the large amount of scatter between data points, 

we cannot rule out the possibility of small values of X /X . As can be seen from 

s u 

Figure 9, an increase of 10 percent in the cross sections would reduce these curves 

by 10 percent but would not alter this conclusion. 

A different conclusion has been reached by the Chicago group (Garcia-Munoz 

et al., 1979). These authors require a PLD of the form exp (-x/X ) with X - 5 g/cm 2 

o o 

2 

of hydrogen (6 g/cm ) of interstellar material) and truncated with a linear rise 

2 

between 0 and 1.3 g/cm (hydrogen) to fit both the B/C and Sub-Fe/Fe ratios 
from 500 MeV/nuc to 1 GeV/nuc. At lower energies they require a lower value 
of X Q . The discrepancy between this result ani our analysis Is due mainly to 
differences in the total (as noted earlier) and spallation cross sections adopted. 
This underlines the necessity for new cross section measurements. Another 
difference is in the shape of the injection spectra used. Garcia-Munoz et al. 

, 1979) have adopted (e + 400 MeV/nuc)“ 2 * 6 as opposed to our (E+ 938 MeV/nuc)“ 2,3 : 







i 


I! 
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This difference is only important at low energies (< 500 MeV/nuc) where 

ionization losses and solar modulation are .important. Their spectrum will clearly 
be too eteep at high energies once energy dependent leakage is taken into account. 

Abundance of Individual Iron Secondaries 

Using our standard model, we have calculated the expected abundances 
of the elements Sc, Ti, V, Cr and Mn relative to Fe as i\ function of energy. 

The results are shown in Figures 10 (a)- (c), where they are compared with 
the available data. As can be seen, the predicted abundance ratios are in 
overall reasonable agreement with the data. In particular, the energy 
dependence of abundance ratlcs from the preliminary data from the HEAO-3 
experiment (Copenhagen-Saclaj Collaboration, 1980) is adequately reproduced 

with the exception of the Mn/Fe ratio which is less steep than that predicted 

54 54 

for Mn stable against P decay. That this is due to the survival of Mn 

above a few tens of GeV/nuc and almost complete decay below ~ 1 GeV/nuc Is shown 

5 A 

in Figure 10(a). We also show our prediction for a mean lifetime of JH Mn 
agalnct p decay of 2 x 10^ years (Cause, 1973) and propagation in a medium 
of mean density 0.33 atoras/cn (Wiedenbeck and Greiner, 1980), This gives 
an acceptable fit to the HEA ( )-3 data although a better fit would result for 
a shorter decay time or a lejs dense propagation region. 

Because of the large uncertainties attached to the measurements df the 
spallation ctoss sections, i: is difficult to go further at the present 
time. Clearly there is interesting structure in these cross sections which 
should be observable in the abundance ratios of Fe secondaries. Indeed, If 
a pronounced maximum occurred in a particular cross section at several 
hundred MeV/nuc, it may be possible to obtain a '’direct measurement" of 
the amount of solar modulation by determining the energy at which the 
abundance of that species relative to Fe was greatest. The displacement 
in energy of this peak relative to* that in the spallation cross sections 
would, after correcting for the effect of ionization losses in the inters 
stellar medium, etc., be equal to the mean energy lost in the heliosphere, 
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Detailed measurements of the relative abundances of Individual Fe-secondary 
Isotopes may be moire suitable for this study. This is because features in 
cross sections of Individual isotopes are expected to be more pronounced 
than in cross sections for the production of a fctoup of isotopes which 
contain the superposition of features at different energies. Cosmic ray 
data with sufficiently good statistics to measure ratios of Individual : 

nuclei to iron as a function of energy to 3 percent accuracy should soon 
be available from the HEAQ-3 satellite experiments. This data should be 
able to determine X to an accuracy of *«> 0.3 g/cm . Before such studies 

B 

can be completed, it will be necessary for new and more accurate measurements 
of the energy dependence of the spallation cross sections to be made. 

Injection Spectra 

As noted earlier, the previous interpretation of the observed energy 

spectra (Ormes and Freier, 1978) in terms of injection spectra of the form 
-2.3 

W * , where W is the total energy per nucleon, needs some revision in order 
to take into account ionization losses in the interstellar medium and sole^ 
modulation. In figure 11, we give a survey of recent observations of the 

-2 3 

energy spectrum of C and Fe. Predictions for the standard model (W 

spectrum) with <f> ■ 0, 200 MV and 400 MV and for injection spectra proportional 
-2 3 

to R ’ (where R is rigidity) with <J> * 200 MV are also shown. For the de- 
gree of solar modulation usually assumed (l.e. 4> * 200 MV), the predicted 

-2.3 

spectra are seen to be too flat at low energies for the W * model and too 
-2 3 

steep for the R * model. It seems that the injection spectra must lie between 
these two forms. 

Summary 

We have calculated the composition of cosmic rays above 100 MeV/nuc 
expected for different pathlength distributions and compared our predictions 
with the data available. Our calculations include the effects of the energy 
dependence of the spallation cross sections, ionization losses, solar modu- 
lation and propagation through the atmosphere down to balloon altitudes. 
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We find that an exponential pathlength distribution is able to account for 
the ratios of secondary- to-prlmary nuclei observed in the cosmic radiation from 

2 

100 MeV/nuc to several hundred GeV/nuc provided the escape length is about 7 g/cm 

* 

below 2 GeV/nuc and proportional to E ^ ^*^at higher energies. Below 500 MeV/nuc 
the cosmic ray composition is mcst sensitive to the shape of the injection spectrum 
and solar modulation. The shape of the PLD is best measured .between 500 MeV/nuc 
and 2 GeV/nuc and, with the advent of the HEAO-3 results will he limited by our 
knowledge of both the partial and total cross sections, particularly for inter- 

t 

actions of iron nuclei. At energies above ~ 2 GeV/nuc the energy dependence of 
the escape length dominates. Using our standard model we predict abundances of 
the individual Fe-secondaries which are consistent with the data available and 
generally lie about 10 percent above the preliminary HEAO-3 data. A 10 percent 
increase in the total cross sections would improve the agreement for most of the 
individual elements, and the general trend of the predicted energy dependence 
seems to be tracked by the data,. An extreme of half the cosmic ray matter being 
around the source \ - 2.6 g/cm seems to be ruled out by these data; however, 
any value of X g ^ 1 g/cm a is possible within the limits of the - data, the cross 
sections, and these calculations. 

When better data on Fe-secondary abundances from HEAO-3 become available 
more refined tests will be possible; however, it will be necessary first to 
determine the energy dependence of the (Be+B)/C ratio and hence the energy 
dependence of the escape length with greater precision. In addition, new 
measurements of the energy dependence of* both the total and the spallation 
cross sections of Fe are vital If the sub-iron cosmic ray data are to be used 
to more .accurately determine the shape of the pathlength distribution. 
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ORIGINAL PAGE YS FIGURE CAPTIONS 

OF POOR QUALITY 

Figure 1 . The Interstellar ratios (Be+B)/C and (21<Z<25)/Fe predicted as 
a function of energy for the standard model (described later) . The effect 
of neglecting ionization losses and the energy dependence of the spallation 
cross sections In propagation calculations is Illustrated for particles outside 
the heliosphere (t}> “ 0). 

Figure 2 . The ratios (Be+B)/C and (21<Z£25) /Fe predicted for the standard* 1 

model outside the heliosphere, above the atmosphere (effect of sdlar modula- 

2 

tlon) and at the depths indicated (g/cm ) in the atmosphere. 

Figure 3 . The (Be+B)/C ratio predicted for an exponential pathlength dlstri- 

2 

bution for the mean escape lengths (g/cm ) indicated. The effects of solar 
modulation with <J> = 200 MV are Included. The experimental data have been 
corrected for propagation through the atmosphere. 

Figure 4 . The variatlor of mean escape length with energy based on our 
interpretation of the (I e+B)/C ratio and I he B/C ratio* The scale at the right 
of the figure gives the mean escape length for the most extreme two-zone model 
(discussed later in the text). 

Figure 5 . A survey of data on the ratio ( 175Z<25) /(26_<Z<28) is compared with 
our prediction for the standard model (so]ar modulation with $ = 200 MV 
is included). The sumbols al'e as follows: G, Garcia-Munoz et al., 1977} 

L, Lezniak and Webber, 1978; J, Juliusson, 1974; 0, Orth et al., 1978; 

Filled circle, Simon et al., 1980. 

Figure 6 . A survey of data on the Fe/O ratio is compared with our prediction 

for the standard model. The effect of solar modulation and propagation 

2 

through 3, 6 and 9 g/cm of atttosphere are indicated. 

Figure 7 . The shape of the pathlength distribution in the tWo-zone model 

for various values of X /A. , keeping A fixed. 

s d e 


Figure 8 . The variation of B/C and (Be+B) /C at high energies with X for 

the two-zone model for various values of X /X. (the numbers attached to 

s e 

the curves). 

Figure 9 . The variation with energy of tbu ratios (Be+B)/C and (21<Z<25)/Fe 
for the standard model and for the two-zone model with X /X. - 0,33. Separate 

S 6 

curves are given for $ - 0 , 200 MV and 400 MV. Symbols for the experimental 
data on the (21<_Z<25)/Fe ratio are as follows: M* Maehl et al., 1977; 

L, Lund et al., 1975; F, Freier et al», 1979; B and 3, Benegas et al., 1975; 

S, Scarlett et al., 1978: Small Diamond, Copenhagen-Saclay Collaboration, 

1980, (preliminary HEAO-3 data); Filled circle, Leznlak and Webber, 1978; Large 
Diamond, Israel et al„, 1979; Rectangle, Garcla-Munoz et al., 1979; Dashed 
Rectangle, Garcla-Munoz et al., 1977. 

Figure 10 . The predicted and observed abundance of individual Fe-secondary 
elements: (a) Mn; (b) Cr and Sc; (c) Ti and V. 

Figure 11 . The measured fluxes of Fe and C are compared with the predictions 
for the standard model with <j> - 0, 200 MV and 400 MV. Also given are 
predictions for an injection spectrum which is a power law in rigidity. 

The curves are normalized such that the flux of C at 10 GeV/nuc is 
3.2xl0" 2 m“ 2 s" 1 sr" 1 (GeV/nuc)“ 1 . 


ORIGINAL PAGE M 

of poor QUALrrif; 


21 


TABLE 1 


Cosmic Ray Cross Sections (mb) 

for the Spallation of 

la C Into Ll» 

Be and B 

E(MeV/nuc) 

100 

251 

630 

1584 

3981 

6 L1 

15.8 

15.8 

15.3 

14.4 

14.4 

7 li 

9.8 

11.8 

13.0 

12.6 

12.6 

7 Be 

11.5 

12,0 

11.7 

10.5 

10, .5 

9 Be 

4.67 

5.44 

5.47 

4.60 

4.60 

1 °Be 

1.28 

2.17 

3.06 

3.04 

3.04 

1 °B ' 

37.5 

32.6 

27.2 

22.4 

22.4 

'’B V ' 

82.1 

61.5 

70.0 

62.9 

61.4 
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TABLE 2 


Cosmic Ray Croas Sections (mb) for the Spallation of Fe into Sub-Fe Elements 


E(MeV/nuc) 

100 

251 

630 

1584 

3981 

Sc 

16.8 

15.6 

29.7 

33.3 

27.7 

Ti 

107.1 

77.2 

118.2 

H5.1 

93.4 

V 

96.6 

53.7 

63,6 

51.6 

42.9 

Cr 

408,2 

183,9 

169.9 

118.9 

104.2 

Mn* 

2/1.7 

144.0 

114.0 

94.2 

92.8 


* 64 Mn assumed stable 
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Predicted Cosmic Ray Composition 


SOURCE INTERSTELLAR ABUNDANCE AT GIVEN ENERGY (MeV/nuc) 




100 

316 

1000 

3162 

10000 

Li 

0 

1.07 

1.29 

1.59 . 

1,46 

0.904 

Be* 

0 

0,618 

0,785 

0.960 

0.843 

0.529 

B 

0 

2.04 

2,21 

2.31 

1,86 

1.20 

C 

3.77 

9.90 

8.81 

8.08 

7.39 

6.32 

N 

0,302 

2.03 

2.06 

2.11 

1.79 

1.32 

0 

4.23 

8.67 

7.90 

7.43 

6.98 

6.23 

F 

0,015 

0.148 

0.174 

0.203 

0,184 

0.133 

m 

0,603 

1.36 

1,34 

1.34 

1.27 

1.08 

Na 

0.026 

0.170 

0.191 

0.213 

0.201 

0.151 

Hg 

0.892 

1.61 

1,56 

1.55 

1.49 

1,33 

A1 

0.113 

0.243 

0.251 

0.277 

0.261 

0.222 

SI 

0,872 

1.28 

1.25 

1,28 

1,24 

1.15 

T 

0 

0.021 

0,031 

0.054 

0,051 

0,035 

S 

0.140 

0.230 

0.254 

0.319 

0.307 

0,257 

Cl 

0.001 

0.024 

0.037 

0,066. 

0.059 

0.041 

Ar 

0.028 

0.089 

0,115 

0.376 

0.158 

0.118 

K 

0.001 

0.044 

0,065 

0.014 

0.097 

0.069 

Ca 

0.060 

0.153 

0,183 

0.249 

0,213 

0.171 

Sc 

0 

0,031 

0,039 

0.057 

0.045 

0,033 

Ti 

0.006 

0.132 

0,151 

0,189 

0.144 

0.111 

V 

0.001 

0.070 

0,073 

0.079 

0.058 

0.046 

Cr 

0.031 

0,212 

0,199 

0.188 

0.147 

0.127 

Mnt 

0.023 

0.120 

0.125 

0,119 

0.110 

0.097 

Fe§ 

1.0 

1.0 

loO 

1.0 

1.0 

1.0 

Co 

0 

0.0008 

0,0011 

0.0014 

0,0011 

0.0010 

N1 

0.042 

0.037 

0,037 

0,038 

0.038 

0.039 


* l0 Be assumed stable 
+ 64 Mn assumed stable 
§ Normalization 


R. J, Protheroe and J. F. Ormei : Code 661, NASA/Goddard Space Fligh Center , 

Greenbelt, MD 20771 


G. M. Comstock: Cosmo Science Associates , 21 Erland Road, Stony Brook, NY 11790 
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